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•  Power  management  optimization  of  fuel  cell  hybrid  vehicles  is  studied. 

•  New  formulation  for  the  optimization  of  real-time  controllers  is  proposed. 

•  Real-time  controllers  can  be  optimized  over  many  driving  cycles  at  once. 

•  Three-parameter  PI  controller  surpasses  a  ten-parameter  fuzzy  controller. 
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Fuel  cell  hybrid  vehicles  offer  a  high-efficiency  and  low-emission  substitute  for  their  internal  combustion 
engine  counterparts.  The  hybridization  significantly  improves  the  fuel  economy  of  the  vehicle;  however, 
exploiting  the  hybridization  requires  a  well-designed  power  management  strategy  that  optimally  shares 
the  power  demand  between  the  power  sources.  This  paper  deals  with  the  optimization  of  power 
management  strategy  of  a  fuel  cell/battery  hybrid  vehicle,  both  off-line  and  in  real-time.  A  new 
formulation  of  the  optimization  problem  for  the  real-time  strategy  is  presented.  The  new  approach  al¬ 
lows  the  optimization  of  the  controller  over  a  set  of  driving  cycles  at  once,  which  improves  the 
robustness  of  the  designed  strategy.  The  real-time  optimization  is  applied  to  two  forms  of  real-time 
controllers:  a  PI  controller  based  on  Pontryagin’s  Minimum  Principle  with  three  parameters  and  a 
fuzzy  controller  with  ten  parameters.  The  results  show  that  the  PI  controller  can  outperform  the  fuzzy 
controller,  even  though  it  has  fewer  parameters.  The  real-time  controllers  are  designed  by  simulation 
and  then  validated  by  experiment. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cell/battery  hybrid  vehicles  are  characterized  by  the  use  of 
two  power  sources:  a  fuel  cell  system  (FCS)  as  the  main  power 
source  and  a  battery  as  an  auxiliary  one.  The  addition  of  the  battery 
leads  to  a  significant  improvement  in  fuel  economy  by  restricting 
the  operation  of  the  FCS  to  high-efficiency  operating  points  and  by 
adding  the  possibility  of  regenerative  braking.  Additionally,  the 
cruising  power  demand  is  assigned  to  the  FCS  whereas  the  battery 
helps  to  meet  the  peak  power  demands  and,  therefore,  the  FCS  can 
be  downsized. 

The  power  management  strategy  defines  the  contribution  of  the 
two  power  sources  in  fulfilling  a  given  power  demand.  For  hybrid 
electric  vehicles  (HEVs),  the  main  objective  of  any  power 
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management  strategy  is  to  minimize  the  fuel  consumption,  while 
limiting  the  variation  of  the  state-of-charge  (SOC)  of  the  battery. 
Most  literature  has  dealt  with  the  power  management  of  the  con¬ 
ventional  hybrid  vehicles  (powered  by  an  internal  combustion 
engine  (ICE)  as  a  main  power  source);  however,  the  design  and 
optimization  approaches  can  be  mostly  adapted  to  the  fuel  cell 
hybrid  vehicles  (FCHVs).  In  comparison  to  conventional  HEVs,  the 
optimization  problem  is  further  simplified  for  FCHVs,  since,  firstly, 
an  on/off  operation  of  FCS  is  not  likely  (which  is  a  normal  operation 
of  ICE  and  it  can  be  optimal  even  during  highway  cruising  1  ])  and, 
secondly,  FCHVs  are  usually  equipped  with  a  single-speed  trans¬ 
mission  and,  therefore,  the  need  for  gear-ratio  optimization  is 
eliminated. 

The  optimal  power  management  strategy  can  be  obtained  when 
the  driving  cycle  is  completely  known  to  get  the  so-called  off-line 
optimum.  Dynamic  Programming  (DP)  [2]  is  still  the  most  common 
off-line  optimization  algorithm  for  HEVs,  for  it  guarantees  the  global 
optimum  without  any  assumptions  regarding  the  convexity  and 
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linearity  of  the  system  models;  however,  DP  has  the  drawback  that  it 
requires  a  long  computation  time.  Convex  optimization  has  been 
recently  applied  to  HEVs  optimization  as  a  significantly  faster  alter¬ 
native  to  DP  if  the  system  can  be  described  by  a  set  of  convex  models 
and,  therefore,  the  accuracy  of  the  optimal  solution  (in  comparison  to 
DP)  depends  on  the  accuracy  of  the  approximations  done  to  relax  the 
non-convex  behavior  of  the  system,  if  any  [3,4  .  Pontryagin’s  Mini¬ 
mum  Principle  (PMP)  [2]  has  been  successfully  applied  to  the  power 
management  optimization  of  HEVs  [5],  where,  in  comparison  to  DP, 
it  is  very  fast;  however,  PMP  gives  only  the  necessary  conditions  that 
must  be  satisfied  by  the  global  optimum  and  not  the  sufficient 
conditions  and,  therefore,  the  PMP  solution  may  be  a  local  optimum. 
However,  if  the  solution  of  PMP  is  unique,  then  it  must  be  a  global 
optimum.  The  uniqueness  of  PMP  solution  is  guaranteed  if  the  bat¬ 
tery  model  parameters  are  constant  independent  of  its  SOC  [5]. 

The  drawback  of  the  off-line  optimal  power  management  strat¬ 
egy  is  that  it  is  highly  dependent  on  the  driving  cycle  and,  therefore, 
the  optimal  strategy  over  one  driving  cycle  may  result  in  a  poor 
performance  over  other  ones  in  real-time.  Hence,  designing  a  near- 
optimal  real-time  power  management  strategy  is  still  an  open  issue 
and  a  matter  of  research.  Equivalent  Consumption  Minimization 
Strategy  (ECMS)  has  been  introduced  as  a  near-optimal  real-time 
strategy  [6-8];  however,  the  excellent  performance  of  ECMS  re¬ 
quires  cycle-dependent  parameters  to  be  pre-programmed  and,  in 
Ref.  [9],  it  has  been  shown  that  ECMS  can  be  regarded  as  a  realization 
of  PMP.  Since  the  problem  associated  with  the  off-line  optimization 
techniques  is  their  cycle-sensitivity,  online  strategy  adaptation 
methods  have  been  devised,  where,  for  instance,  the  adaptation  was 
based  on  driving  pattern  recognition  in  Ref.  [10]  and  on  driving  cycle 
prediction  in  Refs.  [11,12].  Stochastic  Dynamic  Programming  has 
been  also  proposed  to  alleviate  the  cycle-sensitivity  of  the  optimal 
control  law  [13].  Another  approach  to  the  design  of  real-time 
strategy  is  the  static  optimization,  where  a  certain  form  of  the 
real-time  controller  is  predefined  and  optimized  either  analytically 
or  numerically.  As  examples,  the  real-time  controller  took  the  form 
of  an  affine  function  of  the  power  demand  in  Ref.  [14],  the  form  of 
Boolean  rules  in  Ref.  [15  and  the  fuzzy  form  in  Refs.  [16,17].  When 
the  real-time  controller  is  optimized  over  one  driving  cycle,  then  its 
robustness  over  other  driving  cycles  becomes  questionable  [16,17]. 
In  contrast,  when  many  driving  cycles  of  different  nature  are 
engaged  in  the  optimization,  the  robustness  over  a  random  driving 
cycle  would  be  more  guaranteed. 

This  paper  illustrates  the  idea  that  the  real-time  controllers  can 
be  optimized  over  a  set  of  driving  cycles  at  once.  A  comparison  is 
made  between  two  forms  of  real-time  controllers:  a  PI  controller 
based  on  PMP  and  a  fuzzy  controller.  The  simulation  results  are 
validated  on  a  small-scale  test  stand. 

The  paper  is  structured  as  follows.  The  vehicle  model  is  intro¬ 
duced  in  Section  2.  This  model  is  used  to  calculate  the  power  de¬ 
mand  of  the  evaluation  driving  cycles.  In  Section  3,  the  power 
demand  is  scaled  down  to  fit  the  test  stand  and  the  model  of  the 
hybrid  power  source  (fuel  cell  and  battery)  is  presented.  The  off¬ 
line  optimization  using  PMP  is  recalled  in  Section  4.  The  optimi¬ 
zation  problem  of  the  real-time  power  management  strategy  is 
formulated  in  Section  5  and  applied  to  two  types  of  controllers:  a 
PMP-based  PI  controller  and  a  fuzzy  controller.  Experimental  vali¬ 
dation  of  the  simulation  work  is  illustrated  in  Section  6.  The  paper 
is  concluded  in  Section  7. 

2.  Vehicle  model 

A  schematic  of  the  fuel  cell/battery  hybrid  vehicle  is  shown  in 
Fig.  1.  The  battery  is  directly  connected  to  the  bus  whereas  the  fuel 
cell  system  is  connected  via  a  DC/DC  converter.  The  bus  feeds  the 
motor  which,  in  its  turn,  drives  the  vehicle  wheels. 


Fig.  1.  Schematic  of  fuel  cell  hybrid  vehicle. 


The  vehicle  model  aims  to  calculate  the  electric  power  required 
by  the  hybrid  power  source  to  meet  a  given  speed  profile,  v(t).  The 
traction  force,  Ft,  required  at  the  wheels  for  a  given  speed  is  the  sum 
of  four  forces: 

Ft  —  Face  +  +  Fr  +  Fg  (1) 

where  Facc  is  the  force  required  to  accelerate  the  vehicle,  Fa  is  the 
force  required  to  overcome  the  aerodynamic  drag,  Fr  is  the  one  to 
overcome  the  rolling  friction  with  the  road  surface  and  Fg  is  the 
gravity  force  acting  on  the  vehicle  when  driving  on  non-horizontal 
roads  (up  or  downhill).  These  forces  are  given  in  (2). 

Facc=  (m  +  ff)  t 

Fa  =  \p  Q lAfV2 

mg(ju o  +  MiF)cos(a)  v  >  0  ^ 

0  v  =  0 

Fg  =  mg  sin(a) 

The  parameters  of  the  vehicle  model  are  listed  in  Table  l.The 
road  angle,  a ,  is  assumed  0  for  all  the  evaluation  driving  cycles. 
Given  the  traction  force,  the  required  torque  at  the  wheel,  Tw  =  rFt, 
and  the  wheel  rotational  speed,  ww  =  v/r,  are  then  obtained.  The 
transmission  line  is  modeled  by  its  efficiency,  and  gear  ratio,  R. 
The  motor  torque,  Tm i,  and  its  rotational  speed,  wm,  is  then  given  by 

t  _  /  W (F^tr)  Fw  >  0 

ml  “  1  TwWF  Fw<0  (3) 

Wj-ji  -  Rv\/y\j 

In  the  dynamic  state,  the  actual  mechanical  torque,  Tm,  pro¬ 
duced  by  the  motor  is  given  by 

Tm  =  Tml+0m^  (4) 

Given  Tm  and  wm,  the  input  electric  power  can  be  calculated 
using  the  efficiency  map  of  the  motor  and  its  driver,  ijm,  as 
expressed  in  (5).  The  efficiency  map  is  shown  in  Fig.  2. 


Table  1 

Parameters  of  the  vehicle  model. 


Parameter 

Value 

Vehicle  mass,  m  (kg) 

1500 

Wheels  inertia,  0W  (kg  m2) 

4 

Motor  inertia,  0m  (kg  m2) 

0.1 

Wheel  radius,  r  (m) 

0.305 

Gravity  constant,  g  (m  s-2) 

9.8 

Friction  coefficients;  go,  gi  (itT1  s) 

0.007,  0.00012 

Aerodynamic  drag  coefficient,  Cd 

0.3 

Air  density,  p  (kg  m-3) 

1.2 

Vehicle  frontal  area,  Af  (m2) 

1.75 

Gear  ratio,  R 

8.9 

Transmission  efficiency,  r/tr 

0.95 
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Speed  (rpm) 

Fig.  2.  Efficiency  map  of  the  motor  and  its  driver.  The  dotted  line  represents  the  torque 
limits.  Data  taken  from  the  ADVISOR  model  of  a  Westinghouse  75-kW  AC  induction 
motor. 


(  TmWm/ T] m  >  0 

^dem  —  \  VmTmWm  Tm  <  0  (5) 

{  0  Tm  =  0 


3.  Hybrid  power  source 

The  optimization  results  will  be  validated  on  a  test  stand 
composed  of  a  fuel  cell  system,  DC/DC  converter,  10-cell  lithium- 
polymer  battery  and  two-quadrant  load.  The  battery  cells  have 
lithium  nickel  manganese  cobalt  oxide  (LiNiMnCoC^)  as  a  cathode 
material  and  graphite  as  an  anode  material. 

The  test  stand  is  supposed  to  simulate  the  vehicle  on  a  small 
scale.  The  hybrid  power  source  of  the  vehicle  has  a  maximum  fuel 
cell  power  of  50  kW,  maximum/minimum  battery  power  of 
±30  kW,  nominal  bus  voltage  of  320  V  and  battery  capacity  of  7  Ah. 
The  test  stand  is  equipped  with  a  1.2  kW  fuel  cell  system,  which  is 
used  to  calculate  the  scaling  factor  of  0.024  (=1.2  kW/50  kW).  The 
battery  power  limits,  battery  energy  and  the  power  profiles 
calculated  from  the  vehicle  model  are  scaled  down  by  the  same 
factor.  So,  the  battery  power  is  limited  in  the  test  stand  to  ±0.72  kW 
(=  ±30  kW  x  0.024).  The  vehicle  battery  has  an  energy  of  2240  Wh 
(=7  Ah  x  320  V)  and,  therefore,  the  battery  of  the  test  stand  should 
have  an  energy  of  54  Wh  (=2240  Wh  x  0.024).  Assuming  a  37  V 
nominal  battery  voltage  in  the  test  stand,  a  1.5  Ah  battery  capacity 
is  required  in  the  test  stand. 


Actually,  the  battery  used  has  a  capacity  of  30  Ah,  which  is  quite 
large  in  comparison  to  the  required  capacity.  So,  only  1.5  Ah  from 
the  entire  capacity  is  considered  around  60%  SOC,  and  the  reported 
experimental  SOC  corresponds  to  the  required  capacity  and  not  to 
the  actual  one.  This  has  an  implication  on  the  battery  model. 
Although  many  battery  models  have  been  developed  [18],  the 
internal-resistance  model  (Fig.  3)  is  still  the  most  often  used  model 
for  HEV  supervisory  control,  due  to  its  simplicity  for  a  quick 
simulation.  In  this  model,  the  battery  is  represented  by  its  open 
circuit  voltage,  V0c»  in  series  with  its  internal  resistance,  R t>.  Two 
values  of  R b  are  usually  distinguished:  one  when  charging  the 
battery  and  the  other  when  discharging  it.  All  the  model  parame¬ 
ters  vary  with  the  SOC  as  shown  in  Fig.  4.  The  parameters  were 
identified  by  an  FIPPC  test  (Hybrid  Pulse  Power  Characterization 
Profile)  [19].  Since  only  small  part  of  the  capacity  will  be  used,  a 
constant  open  circuit  voltage  of  37.37  V  and  a  constant  internal 
resistance  of  37  mQ  are  used  in  the  simulation.  These  values 
correspond  to  60%  SOC.  It  should  be  noted  that  the  use  of  constant 
battery  model  parameters  has  no  significant  effects  on  the  trans¬ 
ferability  of  our  experimental  results  to  the  real  world  applications, 
since,  on  the  vehicle  level,  the  battery  is  operated  over  a  narrow 
SOC  range  (~45%-75%)  and,  over  such  range,  the  open  circuit 
voltage  and  the  internal  resistance  of  the  battery  do  not  vary 
significantly  with  SOC,  as  illustrated  in  Fig.  4. 

The  hydrogen  consumption  rate,  mh,  was  measured  as  a  func¬ 
tion  of  the  output  power  of  the  DC/DC  converter,  PdC,  and  the 
relation  is  modeled  by  a  quadratic  polynomial  (Fig.  5)  over  the 
range  0—750  W,  which  is  a  sufficient  range  for  all  the  evaluation 
driving  cycles  that  will  be  used.  In  comparison  to  the  linear  inter¬ 
polation  between  the  measurement  data,  using  a  differentiable 
model  results  in  a  smoother  control  of  the  system.  The  combined 
efficiency  of  the  FCS  and  DC/DC  converter  is  defined  as  the  ratio  of 
Pdc  to  the  hydrogen  power,  mh  x  LHV,  where  LHV  (=120  kj  g-1)  is 
the  lower  heating  value  of  hydrogen.  A  comparison  between  the 
measured  efficiency  and  the  model  efficiency  is  illustrated  in  Fig.  6. 

4.  Off-line  optimization 

The  off-line  optimal  power  management  strategy  stands  for  the 
strategy  that  moves  the  system  from  an  initial  state  Xo  at  the 
beginning  of  the  driving  cycle  (t  =  0)  to  a  final  state  Xf  at  the  end  of 
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Fig.  5.  Experimental  hydrogen  consumption  rate  versus  output  power  of  DC/DC  con¬ 
verter  with  a  quadratic  polynomial  fit  to  the  data. 


the  driving  cycle  (t  =  tf)  while  minimizing  the  hydrogen  consumed 
over  the  driving  cycle.  Pontryagin’s  Minimum  Principle  (PMP)  is 
used  to  find  the  optimal  solution  2,5,9,12]. 

The  power  demand,  Pdem,  is  the  sum  of  battery  power,  Pb,  and 
the  power  of  the  DC/DC  converter,  PdC. 

^dem  “  ^dc  +  (6) 

The  state-of-charge  of  the  battery  is  used  as  the  only  state  of  the 
system  (i.e.,  x  =  SOC),  which  is  calculated  according  to  (7)  and 
governed  by  the  dynamic  Equation  (8), 

t 

X(t)  =  X0  -  T  J  /„( T)dT  (7) 

0 


=  it  =  (V-  -  >/V&:-4RbU(t))  (g) 

=  f(x(t),u(t)) 

where  u  is  the  battery  power,  which  is  used  as  a  control  variable.  Jb 
and  Qb  are  the  battery  current  and  battery  capacity,  respectively. 
The  battery  power  and  current  are  considered  positive  for  dis¬ 
charging  and  negative  for  charging. 

The  power  of  DC/DC  converter  is  bounded  between  a  minimum, 
Pdc.min,  and  a  maximum,  Pdc.max-  Likewise,  the  battery  power  is 
bounded  between  Pb.min  and  Pb.max-  These  power  constraints 
impose  a  limitation  on  the  control  variable,  u,  as  follows: 

U(t)  =  Lrnax(Pdem(t)  —  ^d^max^fymin)?  min(^dem(0 

—  ^dc^mim  ^Vmax)  J  (9) 

The  cost  function  to  be  minimized  is  the  total  hydrogen  con¬ 
sumption,  M,  and  it  is  given  by: 


M  =  J  mh(u(t),t)dt  (10) 

0 

The  Hamiltonian  of  the  system  is  defined  as: 

H(x(t),u(t),A(t),  t)  =  rhf(u(t),  t)  -  A(t)/(x(t),  u(t))  (11) 

where  X{t)  is  called  the  co-state,  which  is,  in  general,  a  function  of 
time.  The  Minimum  Principle  states  that  the  optimal  control  vari¬ 
able  minimizes  the  Hamiltonian,  that  is: 

H(x°(t),u°(t),A0(t),t)  <H(x°(t),M°(t),t) 

Vte  [0,  tf] ,  VueU(t) 

where  x°(t)  is  the  optimal  state  trajectory  which  results  from  the 
optimal  control  u°(t)  and  the  initial  state  xo,  and  A°(t)  is  the  optimal 
co-state  trajectory,  whose  dynamics  is  governed  by 

A°(t)  =  f  (x°(f),u0(f),A°(t),t) 

=  -A0(t)|(x°(t),u°(t)) 

The  function /in  (8)  is,  in  general,  a  function  of  x,  since  the  V0c 
and  the  Pb  vary  with  x;  however,  as  stated  in  Section  3,  these  two 
parameters  are  assumed  constant  and,  therefore,  /  is  independent 
ofx  in  our  case.  As  a  result,  according  to  (13),  the  optimal  co-state  is 
constant.  In  this  special  case,  where  the  optimal  co-state  is  con¬ 
stant,  the  solution  of  PMP,  if  exists,  is  unique  and  it  is  a  global  op¬ 
timum  [5]. 

For  a  given  co-state,  the  Hamiltonian  is  minimized  at  each  time 
step  over  the  permissible  control  range  to  find  the  optimal  control. 
The  co-state  plays  the  role  of  a  weighting  factor,  so  that  for  a  high 
co-state  the  charging  of  the  battery  is  favored  and  the  final  SOC 
tends  to  be  high,  and  for  a  low  co-state,  the  final  SOC  tends  to  be 
low.  As  a  result,  the  optimal  co-state  depends  on  the  final  SOC,  Xf, 
which  is  a  predefined  parameter  of  the  optimization  problem.  So, 
finding  the  optimal  co-state  which  leads  to  a  certain  Xf  requires  an 
iterative  procedure  over  the  entire  driving  cycle.  The  iteration 
stops  once  the  required  Xf  is  obtained.  Making  use  of  the  last 
remark,  that  Xf  is  a  monotonically  increasing  function  of  A,  a  binary 
search  algorithm  can  be  used  to  find  the  optimal  co-state  within 
few  iterations.  The  search  algorithm  uses  an  initial  search  space 
Umin,  ^max]-  Theoretically,  X  is  limited  between  zero  and  infinity, 
where  a  zero  co-state  minimizes  the  Hamiltonian  at  the  minimum 
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Fig.  7.  Optimal  PMP  co-state  (top)  and  state  (bottom)  trajectories  over  NEDC  with  and  without  SOC  limitation. 


possible  Pdc  (maximum  possible  Pb)  and  an  infinite  co-state  results 
in  the  maximum  possible  PdC  (minimum  possible  Pb).  In  the  search 
algorithm,  the  initial  Amin  is  set  to  zero  whereas  the  initial  Amax  is 
set  to  4  which  is  more  than  enough  for  our  system.  The  algorithm 
starts  by  testing  Aav,  the  average  of  Amin  and  Amax.  If  the  achieved 
final  SOC  is  lower  than  the  required  one,  Amin  is  increased  to  Aav, 
and  if  the  achieved  SOC  is  higher  than  the  required  one,  Amax  is 
decreased  to  Aav.  In  the  next  iteration,  the  procedure  is  repeated 
with  the  new  search  limits.  At  the  end  of  each  iteration,  the  size  of 
the  search  space  is  halved  and,  within  10  iterations,  the  space  size 
is  reduced  to  less  than  0.1%  of  its  initial  value. 

The  SOC  of  the  battery  is  to  be  limited  between  xmin  =  0.45  and 
Xmax  =  0.75  at  any  time  of  the  driving  cycle.  The  PMP,  as  stated,  does 
not  include  the  SOC  limitation,  in  contrast  to  other  optimization 
algorithms  like  Dynamic  Programming.  There  are  several  ap¬ 
proaches  that  can  be  used  to  incorporate  the  state  constraint  into  the 
formulation  of  PMP  [2,20].  The  approach  used  here  was  introduced 
in  Ref.  [21]  and  can  be  summarized  as  follows:  the  optimal  state 
trajectory  is  calculated  first  for  the  entire  driving  cycle.  If  the  solu¬ 
tion  satisfies  the  state  constraints,  then  the  obtained  solution  is  the 
optimum.  Otherwise,  the  load  profile  is  divided  into  two  sub¬ 
profiles  at  the  time  step  where  the  maximum  violation  of  state 
limits  occurs,  and  the  optimization  is  repeated  for  each  sub-profile. 
If  the  upper  state  limit  is  violated,  then  xmax  is  used  as  a  final  state  for 
the  first  sub-profile  and  as  an  initial  state  for  the  second  one. 
Otherwise,  if  the  lower  state  limit  is  violated,  then  xmin  is  used  as  a 
final  state  for  the  first  sub-profile  and  as  an  initial  state  for  the 
second  one.  The  division  is  repeated  within  each  sub-profile  if  the 
state  constraint  is  violated  again,  until  the  state  constraints  are 
satisfied  for  all  sub-profiles.  Fig.  7  illustrates  the  algorithm  used  to 
consider  the  SOC  limitation  in  PMP.  Here  the  objective  is  to  find  the 
optimal  PMP  solution  over  NEDC  driving  cycle  for  an  initial  SOC,  xo, 
of  0.6  and  a  final  SOC,  Xf,  of  0.65.  At  first,  the  optimal  state  and  co¬ 
state  trajectories  are  calculated  over  the  entire  driving  cycle  be¬ 
tween  t  =  0  and  tf  =  1184  s  without  SOC  limitation.  In  this  case,  a 
violation  of  the  upper  SOC  limit  is  observed  and  the  maximum 
violation  occurs  at  the  time  step  tc  =  807  s.  As  a  result,  the  power 
demand  profile  of  NEDC  is  divided  into  two  sub-profiles, 
Pdemi  =  Pdem  ([0,tc[)  and  Pdem2  =  ^dem  ([WfD-  Over  the  time  span 
[0,rc],  the  optimization  problem  is  re-solved  with  an  initial  state  Xo 


and  a  final  state  xmax  =  0.75,  whereas,  over  the  time  span  [tc,tf],  the 
optimization  problem  is  re-solved  for  an  initial  state  xmax  and  a  final 
state  Xf.  In  this  case,  no  further  sub-division  is  required  since  the  two 
optimal  SOC  sub-trajectories  do  not  violate  the  limitations.  The 
calculated  two  co-state  sub-trajectories  and  the  two  state  sub¬ 
trajectories  are  simply  combined  at  tc  to  form  the  PMP  solution 
with  SOC  limitation.  Comparing  the  unconstrained  and  constrained 
co-state  trajectories  in  Fig.  7,  an  interesting  feature  of  the  co-state  is 
observed  when  the  SOC  limits  are  hit.  The  optimal  co-state  over  the 
time  span  [0,tc  [  is  lower  than  that  of  the  unconstrained  solution 
since  the  achieved  final  state  is  lower,  and  the  optimal  co-state  over 
the  time  span  [t0tf\  is  higher  than  that  of  the  unconstrained  solution 
since  the  initial  state  is  lower.  This  explains  the  jump  in  co-state  at  tc. 

To  sum  up,  the  objective  and  constraints  of  the  off-line  opti¬ 
mization  are  summarized  as  follows: 

tf 

min  J  mh(u(t),t)dt 

m  =mt )) 

x0  =  x(0)  (14) 

xf  =  x(tf) 

*min  —  X  <  *max 
^b,min  —  Pb  —  ^b,max 
^dc,min  —  Pdc  —  ^dc,max 


5.  Real-time  optimization 

In  case  of  off-line  power  management  strategies,  the  driving 
cycle  is  assumed  to  be  known  in  advance  and  an  end  SOC  constraint 
is  assumed.  For  perfect  charge  sustenance,  the  final  SOC  is  required 
to  be  the  same  as  the  initial  one  (i.e.,  Xf = Xo).  In  contrast,  for  real-time 
strategies,  the  future  driving  pattern  is  not  known  and,  therefore, 
the  end  SOC  will  not  necessarily  be  the  same  as  the  initial  one. 

When  designing  a  real-time  power  management  strategy,  this 
strategy  should  be  compared  with  other  strategies,  including  the 
optimal  one,  to  evaluate  it.  Since  different  strategies  may  result  in 
different  end  SOC,  the  comparison  based  only  on  the  fuel 


Objective  : 


Constraints  : 
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consumption  is  not  a  reasonable  approach.  So,  an  algorithm  is 
required  to  quantitatively  compare  the  real-time  strategies. 

The  algorithm  used  here,  is  to  refer  the  fuel  consumption  of  any 
strategy  to  the  same  reference  final  SOC,  Xftref.  The  final  SOC  of  0.6  is 
used  as  a  reference  one.  In  this  case,  we  need  to  define  the  equiv¬ 
alent  fuel  consumption. 

5.1.  Equivalent  fuel  consumption 

The  equivalent  fuel  consumption  is  defined  as  the  fuel  con¬ 
sumption  of  the  strategy  under  interest  if  the  final  state  of  charge 
were  the  reference  one.  The  definition  of  the  equivalent  fuel  con¬ 
sumption  is  illustrated  in  Fig.  8.  The  off-line  optimization  problem 
is  solved  for  various  final  SOC,  Xf,  and  the  resulting  fuel  consump¬ 
tion,  Mopt  =/(x f),  is  shown  in  Fig.  8  for  the  driving  cycle  NEDC,  as  an 
example.  The  resulting  off-line  optimal  curve  is  then  used  to  correct 
the  difference  between  Xf  and  Xf>ref  [22].  If,  for  a  real-time  power 
management  strategy,  a  final  SOC  of  Xf  and  a  fuel  consumption  of  M 
are  achieved  (as  represented  by  the  circle),  a  curve  parallel  to  the 
optimal  consumption  curve  is  constructed  starting  from  Xf  and 
ending  at  Xf>ref  (the  end  is  represented  by  a  triangle).  Then  the 
equivalent  fuel  consumption,  Meq,  is  the  y-coordinate  at  the  end  of 
the  curve.  This  is  expressed  in  (15),  where  M0pt,ref  is  the  optimal 
hydrogen  consumption  for  Xf  =  Xf>ref. 


The  approach  used  here  accounts  for  fuel  consumption  and 
charge  sustenance  as  two  separate  performance  indices.  For  a  given 
driving  cycle,  i,  we  define  two  performance  indices  of  the  real-time 
strategy:  one  for  hydrogen  consumption,  JlHO  and  the  other  for 
charge  sustenance,  Jlsoc,  as  follows: 


M1  -M1 

JUc  =  M'  °Pt’  X  10°(%) 


JsOC  —  \X\  -  xf,ref 


(16) 


The  first  performance  index  measures  how  far  the  hydrogen 
consumption  of  the  real-time  strategy  is  from  the  off-line  optimal 
consumption,  whereas  the  second  one  measures  the  deviation 
from  the  perfect  charge  sustaining  performance.  It  is  important  to 
notice  that  the  definition  of  the  performance  indices  does  not 
distinguish  between  positive  or  negative  (x[  -  xf  ref)  if  the  equiva¬ 
lent  consumption  is  the  same. 

For  a  set  of  driving  cycles,  we  define  the  total  performance 
indices  as  follows: 


JhC  —  JV  ^Jhc 
i 

isoc  =  max/soc 


(17) 


Meq  —  M  —  Mopt  (xf )  +  Mopt  ref  (15) 


5.2.  Optimization  problem  formulation 

Any  real-time  power  management  strategy  has  parameters  that 
can  be  optimized  for  a  given  driving  cycle  in  terms  of  hydrogen 
consumption  and  charge  sustenance.  However,  the  optimized  pa¬ 
rameters  for  a  driving  cycle  may  behave  poorly  for  another,  so  it 
would  be  better  to  incorporate  many  driving  cycles  in  the  optimi¬ 
zation  problem.  In  this  case,  the  probability  of  a  poor  performance 
over  a  random  driving  cycle  decreases.  The  objective  of  the  real¬ 
time  power  management  strategy  is  to  achieve  fuel  consumption 
as  close  as  possible  to  the  off-line  optimal  one  while  keeping  the 
SOC  in  an  acceptable  range. 


Final  SOC 

Fig.  8.  Off-line  optimal  hydrogen  consumption  over  NEDC  versus  final  SOC  for  an 
initial  SOC  of  0.6,  and  an  example  calculation  of  the  equivalent  fuel  consumption. 


where  Jhc  is  the  total  performance  index  for  hydrogen  consumption, 
i.e.,  the  average  of  individual  indices  of  all  driving  cycles,  where  N  is 
the  number  of  driving  cycles  incorporated  in  the  optimization.  Each 
set  of  the  controller  parameters  results  in  a  value  for  the  two  per¬ 
formance  indices  and  the  objective  of  the  optimization  is  to  find  the 
parameters  which  minimize  the  fuel  consumption  while  the  final 
SOC  of  all  driving  cycles  lie  within  a  predefined  range.  The  optimi¬ 
zation  problem  can  be  formulated  as  follows: 

Objective :  min/Hc 

Constraints  :  Vte  [o,  t[]  Vi  (18) 

isoc  <  a 

The  first  constraint  states  that  the  SOC  trajectory  of  all  driving  cy¬ 
cles  must  stay  within  the  allowed  limits.  The  second  constraint 
defines  the  maximum  SOC  deviation  from  the  reference  at  the  end 
of  driving  cycles.  The  optimization  problem  can  be  solved  for 
varying  a.  With  increasing  a,  the  window  of  parameters,  over 
which  Jhc  is  minimized,  expands  and,  therefore,  a  smaller  objective 
value  is  obtained. 

5.3.  Real-time  controllers 

Two  real-time  controllers  are  to  be  optimized.  The  first 
controller  is  based  on  PMP  and  has  a  PI  form  with  three  parameters, 
and  the  second  one  is  a  fuzzy  controller  with  ten  parameters. 

The  form  of  the  PI  controller  is  given  in  (19).  It  controls  the  SOC 
around  the  reference,  xref  =  0.6,  and  outputs  the  co-state  which  is 
then  used  in  real-time  when  minimizing  the  Hamiltonian  defined 
in  (11).  This  type  of  controller  was  used  in  Refs.  [23,24],  as  exam¬ 
ples.  The  controller  has  three  parameters  that  can  be  optimized:  the 
initial  co-state,  Ao;  the  proportional  gain,  /Cp;  and  the  integral  gain, 
K i. 

t 

=  A0  +  I<p  (xref  -  x(t))  +K{  j  (xref  -  x(T))dT  (19) 

0 

The  fuzzy  controller  has  two  inputs  (power  demand  and  battery 
SOC)  and  one  output  (power  of  DC/DC  converter).  Trapezoidal 
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membership  functions  are  used  for  the  inputs  and  the  output,  with 
4  parameters  defining  the  membership  functions  of  each  variable: 
SOC  (xi . .  .X4),  Pdem  (yi  •  •  .3/4)  and  PdC  ( Z\ . .  .Z4),  respectively,  as  shown 
in  Fig.  9.  This  controller  has  a  relatively  low  resolution  with  three 
membership  functions  for  each  variable,  which  are  denoted  as  L 
(Low),  M  (Medium)  and  H  (High).  The  rules  of  the  controller  are 
listed  in  Table  2.  The  logic  behind  these  rules  is  as  follows:  PdC  is  set 
to  high  if  the  SOC  is  low  regardless  of  Pdem  level  (low,  medium  or 
high)  in  order  to  charge  the  battery.  If  the  SOC  is  around  the 
reference  (i.e.,  SOC  is  medium),  PdC  follows  the  level  of  Pdem  in  an 
attempt  to  keep  the  SOC  in  the  medium  range.  If  the  SOC  is  high,  PdC 
is  set  to  a  level  lower  than  that  of  Pdem  when  Pdem  is  medium  or  high 
in  order  to  discharge  the  battery.  The  membership  functions  of  SOC 
are  assumed  to  be  symmetrical  around  the  reference  0.6;  i.e., 
knowing  x\  and  X2  gives  X3  and  X4.  So,  the  parameters  engaged  in 
the  optimization  counts  10:  x\  -  X2,  y\  -  3/4,  Zi  -  Z4.  The  parameters 
of  a  fuzzy  controller  may  be  defined  by  various  means.  In  Ref.  [17], 
for  example,  the  fuzzy  controller  has  the  same  resolution  (i.e.,  three 
membership  functions  for  each  variable)  and  is  defined  by  only  five 
parameters,  whereas  a  higher  resolution  is  used  in  Ref.  [16]  with  22 
parameters.  Generally  speaking,  increasing  the  number  of  param¬ 
eters  increases  the  dimension  of  the  optimization  problem  and, 
therefore,  a  better  controller  performance  is  expected;  however, 
this  comes  at  the  expense  of  increased  optimization  time. 

In  contrast  to  the  PI  controller,  the  optimal  fuzzy  logic  controller 
(with  respect  to  the  optimization  problem  defined  in  (18))  does  not 
necessarily  respect  the  system  power  constraints  (9).  So,  it  is 
possible  for  a  set  of  controller  parameters  that  certain  power  de¬ 
mands  are  not  met  by  the  controller  even  though  the  power  source 
is  physically  capable  of  fulfilling  them.  As  a  workaround,  a  dynamic 
saturation  element  is  added  at  the  output  of  the  controller  to  hardly 
force  the  output  to  be  within  the  limits  allowed  by  the  system.  The 
lower,  Pdc>io,  and  upper,  Pdc.up,  bounds  of  the  saturation  element 
varies  with  power  demand  as  follows: 

*  max(Pdem(0  -  ^b,max^dc,min)  (20) 

^dc,up(0  =  min(^dem(0  —  ^min?  ^dc,max) 

Genetic  algorithm  [25]  was  used  to  optimize  the  two  controllers, 
PI  and  fuzzy,  with  a  fitness  function  and  constraints  as  defined  in 
(18)  with  a  =  0.05  (i.e.,  the  final  SOC  of  any  driving  cycle  must  be 
between  0.55  and  0.65).  Five  driving  cycles  are  engaged  in  the 
optimization,  namely,  NEDC,  UDDS,  NYCC,  Japanese  10-15  and 
HWFET,  whereas  US06  is  used  later  to  test  the  robustness  of  the 
controllers.  The  genetic  algorithm  was  run  for  100  generations  and 
the  population  size  was  chosen  to  be  10  times  the  number  of 
controller  parameters,  i.e.,  30  for  PI  and  100  for  fuzzy.  The  opti¬ 
mization  procedure  proceeds  as  follows:  For  each  set  of  parameters 
(i.e.,  individual)  generated  by  the  genetic  algorithm,  the  corre¬ 
sponding  controller  is  tested  over  each  driving  cycle  separately, 
where  all  driving  cycles  are  initiated  with  an  SOC  of  0.6.  If  the 
controller  satisfies  the  constraints  of  the  problem  (18),  the  corre¬ 
sponding  objective  function,  JHg  is  calculated  according  to 


Table  2 

Rule  list  of  the  fuzzy  logic  controller. 


PdC 

SOC 

L 

M 

H 

Pdem 

L 

H 

L 

L 

M 

H 

M 

L 

H 

H 

H 

M 

Equations  (16),  (17).  Otherwise,  a  high  value  is  assigned  to  the 
corresponding  objective  function  in  order  to  cause  the  algorithm  to 
go  away  from  the  selected  parameters  in  the  subsequent 
generation. 

The  results  of  the  optimization  show  that  the  PI  controller  excels 
the  fuzzy  controller.  The  optimal  Jhc  of  the  PI  controller  is  0.28% 
whereas  it  is  0.74%  for  the  fuzzy  controller.  Those  values  of  Jhc  are 
attained  for  Jsoc  of  0.05  (i.e.,  the  maximum  allowed  Jsoc )•  As  defined 
by  Equations  (16),  (17),  the  optimal  Jhc  represents  the  minimum 
average  deviation  from  the  off-line  optimal  fuel  consumption, 
where  the  average  is  evaluated  over  the  driving  cycles  engaged  in 
the  optimization.  The  optimal  PI  controller  parameters  are  as  fol¬ 
lows:  Ao  =  3.37,  I<p  =  3.10  and  I<[  =  4.70  x  10-5.  A  surface  plot  of  the 
optimal  fuzzy  logic  controller,  together  with  the  output  limitation 
defined  in  (20),  is  shown  in  Fig.  10.  The  message  of  these  findings  is 
the  following:  the  PMP-based  PI  controller  with  only  three  pa¬ 
rameters  can  perform,  on  average,  better  than  a  fuzzy  logic 
controller  with  ten  parameters.  However,  the  difference  between 
the  two  controllers  is  so  small,  that  an  experimental  validation  is 
conducted  for  a  fairer  judgment. 


Fig.  9.  Membership  functions  of  the  fuzzy  controller. 
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Fig.  11.  Photograph  of  the  Hardware-in-the-Loop  test  stand. 


6.  Experimental  validation 

The  results  obtained  by  the  simulation  and  optimization  are 
experimentally  validated  in  this  section,  both  in  open  loop  and 
closed  loop.  The  off-line  optimality  is  validated  in  open  loop, 
whereas  in  closed  loop,  the  performance  of  the  two  optimized  real¬ 
time  controllers,  PI  and  fuzzy,  is  evaluated  on  the  test  stand. 

The  Hardware-in-the-Loop  test  stand  (Fig.  11)  is  composed  of  a 
fuel  cell  system,  DC/DC  converter,  lithium-polymer  battery  and 
two-quadrant  load  to  simulate  both  driving  phase  and  regenerative 
braking.  The  output  current  of  the  DC/DC  converter  is  controlled  by 


Final  SOC 


Fig.  13.  Off-line  optimal  fuel  consumption  over  NEDC  as  a  function  of  final  SOC  in 
simulation  and  experiment. 

the  power  management  strategy.  The  whole  test  stand  is  controlled 
and  monitored  in  real-time  by  an  xPC-Target  computer  equipped 
with  fast  and  high-resolution  input/output  cards.  The  target 
application  is  programmed  in  Simulink  environment  on  a  host 
computer.  The  communication  between  the  target  computer  and 
the  host  computer  serves  to  download  the  compiled  application  to 
the  target  computer  and  to  transfer  the  collected  data  from  the 
target  to  the  host  computer  for  evaluation. 

6.1.  Open  loop  validation 

NEDC  was  used  to  test  the  behavior  of  the  hardware  in  open 
loop.  The  off-line  optimal  co-state  trajectories  calculated  by  the  off¬ 
line  optimization  are  applied  to  the  experimental  power  manage¬ 
ment  strategy  which  includes  the  minimization  of  the  Hamiltonian. 
Four  co-state  trajectories  were  chosen  for  validation  corresponding 
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Fig.  12.  Open  loop  behavior  of  the  test  stand.  (Top)  off-line  optimal  co-state  trajectories  and  (bottom)  the  resulting  state  trajectories  in  simulation  and  experiment.  The  subscripts  1 
and  2  link  the  co-state  trajectories  in  the  upper  sub-figure  with  the  resulting  SOC  trajectories  in  the  lower  one. 


F.  Odeim  et  al.  /  Journal  of  Power  Sources  252  (2014)  333-343 


341 


Table  3 

Experimental  results  of  the  PI  and  fuzzy  controllers  over  the  evaluation  driving  cycles. 


NEDC 

UDDS 

NYCC 

Japan.  10-15 

HWFET 

US06 

PI 

Fuzzy 

PI 

Fuzzy 

PI 

Fuzzy 

PI 

Fuzzy 

PI 

Fuzzy 

PI 

Fuzzy 

Xf 

0.558 

0.568 

0.638 

0.548 

0.649 

0.646 

0.626 

0.579 

0.551 

0.569 

0.608 

0.546 

M(g) 

2.472 

2.578 

3.032 

2.764 

1.046 

1.055 

1.697 

1.548 

3.015 

3.135 

1.631 

1.449 

Meq  (g) 

2.612 

2.686 

2.907 

2.936 

0.888 

0.905 

1.612 

1.617 

3.193 

3.249 

1.603 

1.633 

Fig.  14.  Experimental  co-state  (top)  and  state  (bottom)  trajectories  for  the  PI  controller  over  Japanese  10-15  and  HWFET.  The  simulation  off-line  optima  are  also  shown. 


to  final  SOC  of  0.55,  0.6,  0.65  and  0.7.  Fig.  12  shows  the  first  and 
fourth  co-state  trajectories  together  with  the  corresponding 
simulation  and  experimental  state  trajectories.  The  other  two  tra¬ 
jectories  are  dropped  out  for  the  sake  of  clarity.  As  can  be  seen  from 
Fig.  12,  the  optimal  co-state  makes  a  step-up  when  the  state  tra¬ 
jectory  hits  the  upper  state  limit  of  0.75  (as  discussed  in  Section  4). 


The  small  discrepancy  in  state  trajectory  between  simulation 
and  experiment  results  from  the  instrumentation  error,  mainly, 
from  the  error  in  following  the  load  profile.  Fig.  13  shows  the  cor¬ 
responding  experimental  hydrogen  consumption,  which  is  slightly 
higher  than  the  simulation  optimum  due  to  the  error  in  modeling  of 
the  fuel  cell  system  and  DC/DC  converter  and,  mainly,  from  the 
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Fig.  15.  State  trajectories  of  the  US06  for  the  PI  (top)  and  fuzzy  (bottom)  controllers  with  the  corresponding  off-line  optimum. 
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Fig.  16.  (Top)  US06  power  profile  and  (bottom)  the  output  power  of  the  DC/DC  converter  for  the  PI  and  fuzzy  controllers. 


transient  efficiency  which  is  not  taken  into  account  in  the  models. 
As  will  be  illustrated  in  the  next  subsection,  the  open  loop  evalu¬ 
ation  is  required  if  the  real-time  strategy  is  to  be  compared  with  the 
off-line  optimum. 

6.2.  Closed  loop  validation 

The  two  controllers,  the  PMP-based  PI  and  fuzzy  controllers 
optimized  in  Section  5.3,  were  tested  on  the  test  stand  in  terms  of 
fuel  consumption  and  achieved  final  SOC,  and  the  results  are  listed 
in  Table  3.  The  equivalent  consumption,  Meq,  is  calculated  as 
described  in  Section  5.1.  The  results  show  first  that  both  controllers 
obey  the  design  criterion  that  the  SOC  at  the  end  of  each  driving 
cycle  must  stay  within  ±0.05  from  the  reference  SOC  of  0.6.  The 


Final  SOC  -  Initial  SOC 


Fig.  17.  Experimental  hydrogen  consumption  for  the  off-line  optimum,  the  PI 
controller  and  the  fuzzy  controller  versus  the  difference  between  the  final  and  initial 
SOC,  with  the  corresponding  linear  polynomial  fits,  “opt.”  stands  for  “optimum”  and 
“meas.”  stands  for  “measurement.” 


fuzzy  controller  deviates  a  bit  from  the  criterion  for  UDDS.  The 
second  finding  is  that  the  PI  controller  excels  the  fuzzy  controller  in 
terms  of  fuel  consumption.  The  improvement  in  fuel  consumption 
ranges  from  0.3%  for  Japanese  10-15  to  2.8%  for  NEDC. 

The  idea  behind  the  use  of  a  PI  controller  for  the  co-state  is 
illustrated  in  Fig.  14,  which  shows  the  experimental  co-state  and 
state  trajectories  for  Japanese  10-15  and  HWFET.  The  same 
controller  works  for  these  two  different  driving  cycles,  so  that  the 
SOC  trajectory  respects  the  SOC  limits  and  the  final  SOC  is  in  the 
range  0.55-0.65,  which  are  the  two  design  constraints  of  the 
controller.  If  no  co-state  adaptation  is  used  (i.e.,  the  proportional 
and  integral  parameters  are  set  to  zero  and,  therefore,  the  initial  co¬ 
state  Ao  =  3.37  is  kept  constant  over  the  entire  time  horizon),  the 
simulation  shows  that  HWFET  will  end  up  at  SOC  of  0.23  and  Jap¬ 
anese  10-15  will  end  up  at  SOC  of  0.75,  if  no  SOC  limitations  are 
considered.  This  means  that  a  constant  co-state  of  3.37  is  too  small 
for  HWFET  and  too  high  for  Japanese  10-15;  however,  online 
adaptation  of  the  co-state  via  a  feedback  controller  of  the  SOC 
works  for  different  driving  cycles. 

The  performance  of  the  two  controllers  was,  additionally,  eval¬ 
uated  over  US06,  where  US06  had  not  been  included  in  the  opti¬ 
mization  of  the  controllers.  The  speed  profile  of  US06  was  scaled 
down  to  a  maximum  of  80  km  h-1  to  fit  our  vehicle  model.  The 
results  for  US06  are  reported  in  the  last  two  columns  of  Table  3. 
Also  here  the  PI  controller  saves  about  1.8%  of  hydrogen  in  com¬ 
parison  to  the  fuzzy  controller.  Fig.  15  illustrates  a  comparison  be¬ 
tween  the  off-line  optimal,  the  simulation,  and  the  experimental 
state  trajectories  of  the  US06.  The  simulation  reproduces  the 
experimental  results  very  well  for  both  controllers.  The  PI 
controller  results  in  a  state  trajectory  closer  to  the  off-line  optimum 


Table  4 

Comparison  between  the  pure  fuel  cell  drive  and  the  hybrid  drive  in  terms  of 
hydrogen  consumption. 


Mfc  (g) 

(Mfc  -  Mpi )/MFC  x  100  (%) 

NEDC 

3.267 

20 

UDDS 

3.869 

25 

NYCC 

1.214 

27 

Japan.  10—15 

2.106 

23 

HWFET 

3.636 

12 

US06 

2.106 

24 
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than  the  fuzzy  controller,  which  gives  the  former  the  advantage  in 
terms  of  hydrogen  consumption.  The  power  contribution  of  the  fuel 
cell  side  of  the  hybrid  power  source  is  shown  in  Fig.  16.  The  rest 
power  is  met  by  the  battery. 

A  more  comprehensive  comparison  between  the  PI  and  fuzzy 
controllers  is  illustrated  in  Fig.  17  for  US06.  The  fuel  consumption  of 
the  two  controllers,  together  with  the  off-line  optimum,  is  evalu¬ 
ated  as  a  function  of  the  difference  between  the  final  and  initial 
SOC.  For  the  off-line  optimum,  the  measurements  are  collected  for 
different  final  SOC  between  0.55  and  0.65  in  open  loop  as  described 
in  Section  6.1.  In  closed  loop,  the  initial  SOC  was  varied  between 
0.55  and  0.65  for  the  PI  controller,  and  between  0.5  and  0.6  for  the 
fuzzy  controller.  Linear  polynomial  fits  to  the  data  are  shown  in 
Fig.  17  for  an  easy  comparison.  The  fuel  consumption  of  the  PI 
controller  is  barely  distinguishable  from  the  off-line  optimum.  The 
intercept  of  the  fit  line  represents  the  hydrogen  consumption  for 
the  perfect  charge  sustenance  (i.e.,  final  SOC  =  initial  SOC).  The 
charge  sustenance  is  achieved  by  the  PI  controller  for  an  initial  SOC 
of  0.608  and  the  corresponding  hydrogen  consumption  is  1.616  g, 
which  is  exactly  the  same  consumption  of  the  off-line  optimum. 
The  fuzzy  controller  achieves  the  perfect  charge  sustenance  for  an 
initial  SOC  of  0.546  with  a  consumption  of  1.639  g.  So,  the  PI 
controller  saves  about  1.4%  of  hydrogen  in  comparison  to  the  fuzzy 
controller  when  both  controllers  achieve  the  charge  sustenance. 
More  accurately,  for  a  quadratic  polynomial  fit  to  the  data,  the 
calculated  hydrogen  saving  becomes  1.7%,  which  is  almost  the  same 
as  the  one  reported  in  Table  3.  Lastly,  in  comparison  to  the  single¬ 
valued  evaluation  used  in  Table  3,  the  experimental  evaluation 
illustrated  in  Fig.  17  gives  a  more  comprehensive  overview  of  the 
performance  of  the  real-time  controllers;  however,  it  requires  quite 
extensive  experimental  measurements. 

The  advantage  of  hybridization  with  respect  to  hydrogen  con¬ 
sumption  was  evaluated  by  measuring  the  hydrogen  consumption 
when  only  the  fuel  cell  system  is  used  to  drive  the  load;  however, 
the  battery  was  kept  connected  to  the  bus  to  absorb  the  recuper¬ 
ated  energy.  The  measured  consumptions,  MFc,  over  the  evaluation 
driving  cycles  are  listed  in  Table  4.  Table  4  also  lists  the  percentage 
hydrogen  saving  when  the  hybrid  power  source  is  controlled  by  the 
optimized  PI  controller  in  comparison  to  the  pure  fuel  cell  drive, 
where  the  hydrogen  consumption  of  the  PI  controller,  Mpi,  is  the 
equivalent  consumption  taken  from  Table  3.  Over  most  of  the 
evaluation  driving  cycles,  the  hydrogen  saving  ranges  between  20% 
and  27%.  The  hybridization  is  less  beneficial  for  HWFET,  with  only 
12%  hydrogen  saving,  due  to  the  short  regenerative  braking  phases 
for  this  highway  cycle.  Lastly,  the  values  given  in  Table  4  should  be 
cautiously  considered,  because  several  factors,  such  as  the  change 
in  fuel  cell  system  size  and  vehicle  mass,  are  not  taken  into  account 
and  they  may  affect  the  hydrogen  saving  of  the  hybrid  vehicle  as 
compared  to  the  non-hybrid. 

7.  Conclusion 

The  optimization  of  the  power  management  of  a  fuel  cell  hybrid 
vehicle  was  investigated  through  simulation  and  experimentation. 
The  off-line  optimization  using  Pontryagin’s  Minimum  Principle 
was,  first,  recalled  to  calculate  the  best  hydrogen  consumption  that 
can  be  achieved  when  the  complete  load  profile  is  known  in 
advance.  The  off-line  optimum  was  then  used  as  a  benchmark 
when  optimizing  the  real-time  power  management  strategy,  where 
both  hydrogen  consumption  and  charge  sustenance  were  taken 
into  account  as  two  separate  performance  indices.  It  was  shown 
that  the  real-time  strategy  can  be  optimized  over  a  set  of  driving 
cycles  and  not  only  one,  which  increases  the  robustness  of  the 
designed  controllers  when  they  are  applied  to  real  driving 


conditions.  Two  real-time  controllers  were  tested:  a  PMP-based  PI 
controller  with  three  parameters  and  a  fuzzy  logic  controller  with 
ten  parameters.  Both  simulation  and  experiment  showed  that  the 
PI  controller  can  achieve  a  better  performance  in  terms  of  hydrogen 
consumption,  even  though  it  has  fewer  parameters  to  be  optimized. 

The  performance  of  the  optimized  real-time  controller  depends 
on  two  main  factors:  the  form  of  the  controller  and  the  number  of 
parameters  engaged  in  the  optimization.  For  more  parameters,  a 
better  performance  would  be  expected  when  the  chosen  parame¬ 
ters  are  independent.  Increasing  the  number  of  membership 
functions  of  the  fuzzy  controller  would  lead  to  a  better  performance 
in  terms  of  hydrogen  consumption  and  robustness.  Flowever,  firstly, 
achieving  a  significant  improvement  in  comparison  to  the  simple  PI 
controller  is  doubtful  and,  secondly,  increasing  the  number  of 
optimization  parameters  would,  necessarily,  require  longer 
computation  time  in  order  to  guarantee  the  global  optimum. 
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